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ABSTRACT 

Although they are only minor constituents of the interstellar medium, 
halogen-containing molecules are of special interest because of their unique ther- 
mochemistry. Here, we present a theoretical study of the chemistry of interstellar 
molecules containing the halogen elements chlorine and fluorine. We have mod- 
eled both diffuse and dense molecular clouds, making use of updated estimates for 
the rates of several key chemical processes. We present predictions for the abun- 
dances of the three halogen molecules that have been detected to date in the inter- 
stellar medium: HF, CF+ and HCl. As in our previous study of fluorine-bearing 
interstellar molecules, we predict HF to be the dominant gas-phase reservoir of 
fluorine within both diffuse and dense molecular clouds; we expect the Herschel 
Space Observatory to detect widespread absorption in the HF J = 1 — transi- 
tion. Our updated model now overpredicts the CF"*" abundance by a factor ^ 10 
relative to observations of the Orion Bar; this discrepancy has widened because 
we now adopt a laboratory measurement of the CF"*" dissociative recombination 
rate that is smaller than the estimate we adopted previously. This disagreement 
suggests that the reaction of C"^ with HF proceeds more slowly than the cap- 
ture rate assumed in our model; a laboratory measurement of this reaction rate 
would be very desirable. Our model predicts diffuse cloud HCl abundances that 
are similar to those predicted previously and detected tentatively toward ( Oph. 
Two additional species are potentially detectable from photodissociation regions: 
the H2CI"'" and HCl"*" molecular ions. Ortho-H2Cl"'" has its lowest-lying transi- 
tion in the millimeter spectral region observable from the ground, and the lowest 
rotational transition of HCl"'" is observable with HerscheFs HIFI instrument. 

Subject headings: ISM: Molecules — ISM: Abundances — ISM: Clouds - molec- 
ular processes - infrared: ISM - submillimeter 
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1. Introduction: the unique thermochemistry of fluorine and chlorine 

Of the ~ 150 distinct molecules detected thus far in interstellar clouds and circumstellar 
outflows, over 90% contain only hydrogen and/or elements in groups IVA, VA, and VIaEI 
of the periodic table: viz. carbon, nitrogen, oxygen, silicon, phosphorus, and sulfur. Of the 
remaining 12 detected molecules, eight contain elements in groups lA, llA and lllA (viz. 
sodium, potassium, magnesium and aluminum) and are found exclusively in circumstellar 
envelopes (e.g. Ziurys 1996, and references therein) with abundances that are set in cool 
stellar photospheres; one contains iron (FeO, detected toward the Sgr B2 interstellar gas 
cloud by Walmsley et al. 2002); and the remaining 3 are interstellar molecules containing 
elements of the halogen group (VllA): HCl, HF, and CF"*". 

Although they are only minor constituents of the interstellar medium, halogen-containing 
molecules are of special interest because of their unique thermochemistry. The basic behav- 
ior of any atom in the interstellar medium is determined by two or three key thermochemical 
considerations. First, the flrst ionization potential (IP) determines whether the atom will 
be neutral or ionized within diffuse molecular clouds: an atom X with IP(X) > IP(H) = 
13.6 eV will be shielded by atomic hydrogen from the interstellar radiation fleld (ISRF) and 
will be predominantly neutral, while an atom with IP(X) < 1P(H) will be predominantly 
ionized. Second, the dissociation energy of the neutral hydride, Do(HX), determines whether 
X can react exothermically with H2: atoms for which Dq{HX.) > Do{B.2) = 4.48 eV can react 
exothermically with H2 to form HX, while atoms with Do(HX) < Z}o(H2) cannot. Third (and 
of importance mainly for atoms with IP < 13.6 eV), the dissociation energy of the hydride 
molecular ion, Do(HX+), similarly determines whether the ion X"*" can react exothermically 
with H2 to form HX"*". Fluorine and chlorin^ are unique in regard to the second and third 
criteria, each in a different way: F is the only atom - and CI"*" is the only ion with an ap- 
pearance potential < 13.6 eV - that can react exothermically with H2 to form a diatomic 
hydride. 

The uniqueness of these elements are presented graphically in Figure 1, which shows 
DoiHX.) for each atom X in the flrst three rows of the periodic table. The molecular binding 
energy of HX is largest for fluorine and decreases, nearly monotonically, as one moves left- 



^Hcre we use the traditional designations used in the United States. The current lUPAC designations for 
groups IVA, VA, and VIA are respectively 14, 15, and 16; those for groups lA, IIA, IIIA and VIIA are 1, 2, 
13, and 17 

^In this paper we combine our attention to CI and F; bromine and iodine have much smaller cosmic 
abundances than chlorine and fluorine, and molecules containing bromine and iodine are unlikely to be 
detectable with current technology. 
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wards and downwards in the periodic table. In Figure 2, Do(HX) and -Do(HX"'") are shown 
for all the elements of which interstellar molecules have been detected. Red squares apply 
to atoms with IP < IP(H) and blue squares to those with IP > IP(II). With the exception 
of the halogen elements, all the elements fall into two classes: 

(1) Oxygen and nitrogen, with IP > IP(H), Dq(HX) < Dq(H2) andDo(HX+) > Dq(H2). 
These elements are predominantly neutral in the diffuse ISM. Their atoms do not react with 
H2 except at elevated temperatures that can drive endothermic reactions, and their chemistry 
is driven by cosmic ray ionization. 

(2) Carbon, silicon, sulfur, and phosphorus, with IP < IP(H), Dq(HX) < Dq(H2) and 
Dq(HX~^ ) < Dq(H2). These elements are predominantly ionized in the diffuse ISM. Their ions 
do not react with H2 except at elevated temperatures that can drive endothermic reactions. 

Chlorine and fluorine form two additional classes with one member each: 

(3) Fluorine, with IP > IP(II) and Dq(HX) > Dq(H2). Fluorine is predominantly 
neutral in the diffuse ISM, and can react exothermically with H2 to form HF. 

(4) Chlorine, with IP < IP(H), Dq(HX) < Do(H2) and Do(HX+) > Dq(H2). Chlorine 
is predominantly ionized in atomic clouds. Cl^ can react exothermically with H2 to form 
HC1+. 

The chemistry of interstellar fluorine has been considered in a previous study by Neufeld, 
Wolfire & Schilke (2005; hereafter NWS05), while that of chlorine has been discussed in 
several studies, most recently those of Schilke, Wang & Phillips (1995) and Amin (1996). In 
the present paper, we extend and update these previous investigations to account for recent 
developments in experimental and theoretical studies of relevant chemical processes. Our 
study is also motivated by the prospects for upcoming observations of halogen-containing 
molecules, particularly with the Herschel Space Observatory. 

2. Model 

2.1. Chemical network for chlorine-bearing species 

The key chemical pathways leading to the formation and destruction of chlorine-bearing 
molecules have been elucidated in the papers of Jura (1974), Dalgarno et al. (1974), van 
Dishoeck & Black (1986; hereafter vDB86), Blake, Anicich & Huntress (1986; hereafter 
BAH86), Schilke, Phillips & Wang (1995), Federman et al. (1995), and Amin (1996). In 
diffuse molecular clouds, the chemical network is initiated by reaction of CI"*" with H2: 
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C1+ + H2 ^ Ha+ + H, (1) 

which is followed by dissociative recombination of HCl"^ or a further H atom abstraction 
reaction to form H2C1^: 

HC1+ + H2 ^ H2CI+ + H (2) 
The H2CI"'" molecular ion undergoes dissociative recombination to form either HCl or CI 

H2a+ + e HCl + H (3) 

H2CI+ + e CI + products (4) 

The branching ratio for dissociative recombination of H2CI"'' has not yet been measured, but 
has been assumed in previous studies to be only 10% for production of HCl (reaction 3). 
This assumption was largely motivated by the small HCl abundances observed in diffuse 
clouds (see §4 below.) HCl is destroyed by photodissociation, by photoionization and by 
reactions with He"*", H;^ and C"*". 

In regions where H2 is present, reactions (1), (2) and (4) can rapidly convert Cl^ to 
CI. The rate at which the reaction network is initiated (reaction 1) is therefore limited 
by the availability of CI"*" ions. These are produced primarily by the photoionization of 
CI, following the absorption of ultraviolet photons in the 912 — 958A wavelength range. 
Near cloud surfaces, the UV radiation field is dominated by starlight (i.e. by the interstellar 
radiation field (ISRF) or nearby OB associations), whereas in the deep interiors of dense 
clouds, the radiation largely comprises H2 Lyman and Werner band emissions excited by 
the secondary electrons produced by cosmic rays. Here, proton transfer from to CI, 
producing HCl"*", can also initiate the formation of Cl-bearing molecules. 

Taking as our starting point the reaction list and compilation of rate coefficients in 
the UDFA (formerly UMIST) rate file (Woodall et al. 2007), we have reviewed the rate 
coefficients critically and implemented the following changes: 

1) For reaction of HCl with He"*", H3" and C"*", we have used the statistical adiabatic 
capture model (SACM; Troe 1985, 1987, 1996; see NWS05 for details of exactly how this 
was implemented) to estimate the temperature dependence of the capture rate. Over the 
temperature range 10 - 300 K, we found the capture rate to be proportional to 
for rotationally-cold HCl; accordingly we adopted a dependence for the assumed 

reaction rate. Similar considerations lead to a 7"^°'^^^ temperature dependence for the rate 
coefficient for the reaction of H2O with H2CI"'". We also adjusted the T = 300 K rate 
coefficients in the UDFA compilation, which were based upon measurements of the thermal 
rate coefficients, to account for the fact that HCl is primarily in J = (i.e. rotationally-cold) 
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under the astrophysical conditions of typical interest; here again, we used the SACM model 
to estimate the necessary correction factors. 

2) As in NWS05, we adopted a rate coefficient of 2 x 10~^(T/300 K)-^/^ ^m^ s"^ for 
the dissociative recombination of diatomic molecular ions (in the absence of experimentally- 
determined values). For the dissociative recombination of H2CI"'", we adopted a rate co- 
efficient of 1.2 X 10~^(T/300 K)"*^'^^ cm^ s~^, based upon preliminary results from a recent 
laboratory measurement (Geppert & Hamberg 2009). Because the branching ratios for the 
production of HCl and CI were not determined in this experiment, we followed previous au- 
thors in adopting - as our standard assumption - a branching ratio of 0.1 for the production 
of HCl. As discussed in §4.2 below, we also investigated how the model predictions depend 
upon the assumed value of this branching ratio. 

3) We adopted the rate coefficients of Pradhan & Dalgarno (1994) for charge transfer 
from H"*" to CI, and for the reverse reaction (from CI"*" to H). These reactions were not 
included in the UDFA compilation, but had been considered by BAH86, who assumed a 
large rate coefficient of 10~^ cm^s~^ for the reaction of H"*" with CI and concluded that this 
process was among the subset of most significant reactions. The calculations of Pradhan 
& Dalgarno (1994), however, imply rate coefficients that are smaller than those guessed by 
BAH86 by factors of ~ 20 and ~ 250 at temperatures of 300 K and 100 K respectively. 

4) We have also included the reaction of CI with H2 to form HCl and H, which is 
slightly endothermic (by ~ 0.05 eV), along with its reverse reaction. This reaction, which 
is not included in the UDFA compilation, has been the subject of intensive investigation, 
both experimental (e.g. Kumaran, Lim & Michael 1994 ; hereafter KLM94) and theoretical 
(e.g. Manthe, Capecchi & Werner 2004 ). Here we adopted the rate coefficients recom- 
mended by KLM94 (their equations 10, 11 and 6), which were obtained by a fit to their own 
measurements together with those obtained in previous studies. 

5) In computing revised estimates of the HCl photodissociation and photoionization 
rates, we have used the photoabsorption cross-sections obtained by Brion, Dyke & Cooper 
(2005) by means of dipole (e, e) spectroscopy using an electron energy loss spectrometer 
(EELS). This technique is believed to offer significant advantages over direct UV absorption 
measurements for the case of strong absorption bands in which line saturation is important 
(Chan, Cooper & Brion 1991) For example, in the case of excitation to the M^H state by 
photons near 110.3 nm, Brion et al. obtained an oscillator strength that was a factor of 
4 larger than that obtained in the earlier UV absorption experiments of Nee, Suto & Lee 
(1986). The continuous absorption cross-sections obtained by Brion et al. (2005) are in good 
agreement with previous UV absorption measurements (e.g. Cheng et al. 2002). 
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To determine the photodissociation rate implied by the measured photoabsorption cross- 
sections, the photodissociation probabihty must also be known. For radiation longward of 
115 nm, the fluorescence yield has been measured (Nee et al. 1986) and found to be small, so 
the photodissociation probability is very close to unity. For absorption in the 105.5 - 115 nm 
region, the fluorescence at wavelengths longer than 115 nm was observed to be weak, but the 
experiment of Nee et al. (1986) did not exclude the possibility of strong fluorescent emission 
in the 105.5 - 115 nm region. For the 97.3 - 105.5 nm spectral region, we are unaware of any 
experimental constraints on the fluorescent yield. In computing the photodissociation rate, 
we have assumed that the dissociation probability is unit} § for all radiation longward of the 
ionization threshold at ~ 97.3 nm. This leads to a photodissociation rate of 1.7 x 10~^s~^ 
for an unshielded molecule exposed to the mean interstellar radiation field given by Draine 
(1978), a value that is ~ 70% higher than that given by van Dishoeck, van Hemert & 
Dalgarno (1982; hereafter vDvHD). This difference is largely accounted for by the fact that 
the total oscillator strength shortward of 129 nm is now known to exceed significantly the 
simple estimate given by vDvHD. We obtain an HCl photoionization rate of 1.0 x 10~^°s~^. 

In considering the effect of shielding upon the HCl photodissociation rate, we used a 
method similar to that adopted by NWS05. Here, we modeled the case of a semi-infinite 
slab that is illuminated isotropically, using the photodissociation region (PDR) model of Le 
Bourlot et al. 1993 and Le Petit et al. 2006 to determine how the UV radiation field varies 
with position. We assumed the photodissociation rate to diminish as E2{kAv), where Ay is 
the visual extinction in magnitudes behind the slab surface, E2 is the exponential integral 
of order 2, and k is an adjustable parameter that we varied in order to obtain the best fit 
to the predicted photodissociation rate. The UV dust albedo assumed here was 0.32, the 
forward scattering function was 0.73 (Li & Draine 2001), and the best fit was obtained for 
k = 2.1. 

6) Our estimate of the CI photoionization rate was based upon measurements of the 
ionization cross-section by Ruscic & Berkowitz (1983 ), which are in good agreement with 
theoretical calculations of Brown, Carter & Kelly (1980 ). Two autoionizing resonances are 
present at wavelengths longer than the Lyman limit, but make a negligible contribution to 
the overall photoionization rate. We derived a photoionization rate of 4.8 x 10~^^s~^ for an 
unshielded chlorine atom exposed to the mean interstellar radiation field given by Draine 
(1978), a value that lies a factor of ~ 2 below that adopted by van Dishoeck (1988). 



■^At the other extreme, if we assume the photodissociation yield to be zero for radiation shortward of 115 
nm, the unshielded photodissociation rate is 0.77 x 10^^ s^^. This therefore constitutes a lower limit on the 
true value, and is in good agreement with the corresponding minimum value given by vDvHD. 
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Because the IP of atomic chlorine is only slightly less than that of H, there is only a 
narrow wavelength range, 912 - 958A, over which photoionization can occur. Radiation in 
this spectral region is subject not only to absorption by dust, but also to absorption by the 
H2 Lyman and Werner bands and the HI Lyman lines. While the effect of HI absorption is 
small, that of H2 absorption can be significant. We have considered dust and H2 absorption 
independently, by running the Le Petit et al. (2006) PDR model both with and without 
H2 absorption. We have found that the reduction in the CI photoionization rate can be 
represented to good approximation by a factor 

exp(-iV(H2)/4.28 x lO^^cm-^) 

E2{6.bAv) X , (5) 

1 + v/iV(H2)/4.9 X 1020cm-2 

where A''(H2) is the column density of H2 to the cloud surface. 

The photoionization of HCl occurs over a similarly narrow range of wavelengths. For 
that process, an identical treatment yields a reduction factor 

E.(3A5 A.) X exp(-iV(H,)/3.5 x 

(T/300K)1 + v/iV(H2)/3.2 x 1020cm-2 

7) Using the cross-sections cited in (5) and (6) above, we have computed the rates 
of HCl photodissociation, HCl photoionization, and CI photoionization resulting from H2 
Lyman and Werner band emissions excited by secondary electrons produced by cosmic rays. 
Here we used the same methodology adopted in NWS05. 

8) For the radiative recombination of CI"*", we adopt the value recommended by Mazzotta 
et al. (1998). Dielectronic recombination is negligible at the temperatures of present interest. 

Our model does not include the photodissociation of molecular ions containing the 
halogen elements. With the exception of HCl^, for which photodissociation cross-sections 
have been presented by Pradhan, Kirby & Dalgarno (1991) and by Andric et al. (2002), we 
are not aware of any estimates of the relevant cross-sections. The HCl"*" photodissociation 
cross-sections imply an unshielded photodissociation rate of 5.6 x 10~^^xuv^~^ ■ In the region 
where HCl"*" is most abundant, photodissociation makes a negligible contribution - relative 
to dissociative recombination - to the HC1+ destruction rate; we assume photodissociation 
to be similarly negligible for CF"*" and H2CI"'". 

The full reaction list for chlorine-bearing species is given in Table 1. It contains 21 
processes involving the following species: CI, HCl, CI"*", HCl"*", H2CI"'" and CCl"*". 
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2.2. Chemical network for fluorine-bearing species 

The chemical network for fluorine-bearing species has been discussed by NWS05. The 
chemistry is initiated by reaction of H2 with F: 

F + H2 ^ HF + H (7) 

Once the H2 abundance becomes appreciable, HF becomes the dominant reservoir of fluorine; 
it is destroyed slowly by photodissociation, and by reactions with He"*", H3" and C"*". Reaction 
with C"*" leads to the formation of CF"*": 

HF + C+ ^ CF+ + H (8) 

The fluoromethylidynium ion, CF"*", is destroyed rapidly by dissociative recombination, but 
can nevertheless account for as much as ~ 1% of the gas-phase fluorine abundance. It has 
recently been detected in the interstellar medium (Neufeld et al. 2006). 

For completeness, we have implemented flve updates to the chemical network presented 
by NWS05. Of the following changes, the flrst three of which were noted in the UDFA 
reaction list (Woodall et al. 2007), only the flrst is expected to have any signiflcant effect on 
the model predictions: 

1) For the rate coefficient for dissociative recombination of CF+, we now adopt the recent 
laboratory measurement of Novotny et al. (2005). The adopted rate coefficient lies factors 
of 4 and 1.5 below that assumed by NWS05 at temperatures of 300 and 10 K respectively. 

2) We have corrected an erroneous value adopted by NWS05 for the rate coefficient for 
reaction of F with H2O. The correct value is 1.6 x 10~^^ cm~^s~^. This was stated correctly 
in the text of NWS05, but a value 10 times too large was erroneously given in the reaction 
list table (and erroneously implemented in our calculation). 

3) For the reaction of H2 with F"*", we include a channel that produces Hj' in addition 
to that which produces HF"*". 

4) For the cosmic ray induced photodissociation of HF, we include the effects of H 
Lyman alpha radiation. This results in a 50% increase in the adopted photodissociation 
rate. 

5) For the radiative recombination and dielectronic recombination of F"*", we adopt the 
rate coefficients of Badnell (2006) and Badnell et al. (2003). 

The full reaction list for ffuorine-b earing species is given in Table 2. It contains 18 
processes involving the following species: F, HF, F"*", HF"*", H2F"'", CF"*" and SiF"*". 
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2.3. Interstellar cloud models 

As in NWS05, we have used a modified version of the PDR models of Kaufman et al. 
(2006), Kaufman et al. (1999), Wolfire et al. (1990), and Tielens & Hollenbach (1985) to 
investigate the chemistry within diffuse and dense molecular clouds. The modifications to 
these models have been largely described in NWS05, with the exception of the following: 

1) For the radiative and dielectronic recombination of C^, we adopt the recently- 
computed rate coefficients of Badnell et al. (2006) and Badnell et al. (2003). 

2) In treating the depletion of fluorine and chlorine that results from freeze-out onto 
grain mantles within dense molecular clouds, we use equation (5) of NWS05, but with the 
HF or HCl photodissociation rate augmented so as to include the effects of H2 Lyman and 
Werner band emissions excited by the secondary electrons produced by cosmic rays. As in 
NWS05, we assume a gas-phase fluorine abundance of 1.8 x 10~^ relative to H nuclei in diffuse 
molecular clouds. For chlorine, we assume an analogous abundance of 1.8 x 10~^, the average 
value observed in three sightlines that intersect diffuse molecular clouds: ( Oph (Federman et 
al. 1995), HD 192639 (Sonnentrucker et al. 2002) and HD 185418 (Sonnentrucker, Friedman 
& York 2006). 

In our standard models, as in NWS05, we assume a primary cosmic-ray ionization rate 
of (p = 1.8 X 10"^^ s~^ per H nucleus. Recent observations of Hj]" have led to the postulate 
of substantially larger rates in diffuse clouds (e.g. Indriolo et al. 2007). To investigate the 
effect of enhanced cosmic ray ionization rates upon the chemistry of interstellar halogens, 
we have also considered models for which the assumed (p is increased by a factor 10. 

Our standard models also compute the thermal balance and equilibrium gas temperature 
as a function of distance into the cloud. To investigate the possible effects of enhanced 
temperatures (such as might result from the passage of interstellar shocks), we have also 
considered a series of models in which the temperature is set to some depth-independent 
value (larger than that computed by our self-consistent treatment of the thermal balance). 

For our standard models, we have followed previous authors in adopting a branching 
ratio of 0.1 for the production of HCl following dissociative recombination of H2CI"'". To study 
the dependence upon this branching ratio, we have also considered diffuse cloud models in 
which values of 0, 0.3 and 1.0 were adopted. 
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3. Results 

We have constructed a series of one-sided models in which ultraviolet radiation is in- 
cident upon a semi-infinite slab, with the specific intensity constant over the incoming 
hemisphere. Each model is characterized by two key parameters: the density of H nu- 
clei, TT-H = n{}i) + 2 71(112), assumed constant throught the slab; and the incident radiation 
field, Xuv, expressed in units of the mean ISRF given by Draine (1978). For one-sided slab 
models, we have considered all combinations of nn = 10^, 10^, 10^, 10^, and 10''cm~^ and 
Xuv = 10^ 10^, 10^ and 10^ 

In Figure 3, we present example results for a single such combination: rin = lO^cm"^, 
and Xuv = 10^- Here, we plot the relative abundances of H2, CI, CI"*", HCl, HCl^, H2C1^, HF, 
and CF^ as a function of the depth into the cloud. CCl"^ fails to appear on this plot because 
its abundance is too low. Each abundance is normalized relative to the relevant gas-phase 
elemental abundance (i.e. hydrogen for the case of H, chlorine for the case of Cl-containing 
molecules, fiuorine for the case of F-containing molecules). The depth is measured in terms 
of the visual extinction. Ay, in magnitudes, along a perpendicular ray to the cloud surface. 
Figure 4 shows the gas temperature and the F and CI gas-phase elemental abundances as a 
function of depth. 

In Figures 3 and 4, solid lines show the results for an assumed primary cosmic ray 
ionization rate of 1.8 x 10~^^s~^ per H nucleus (standard case), while dashed lines show the 
results for an assumed primary cosmic ray ionization rate of 1.8 x 10~^^s~^ per H nucleus. 
The effect of the increased cosmic ray ionization rate is negligible for Ay < 3, because the 
ionization of CI is dominated by the ISRF near the cloud surface. At larger depths below 
the cloud surface, the predicted abundances of H2CI"'" and CF"*" are significantly enhanced 
if a value of 1.8 x 10^^^ s~^ is adopted for (p] however, at these depths, the observed H3 
abundances do not support an enhanced value for the cosmic ray ionization rate (Indriolo et 
al. 2007). 

By integrating the molecular abundances into a depth Ay = 10 mag, we have obtained 
the total column densities for each case. In Figure 5, we present the column densities of 
HCl, HCl"*", H2CI"'", HF, and CF"*". Here, we have plotted the molecular column densities as 
a function of nn. Results are shown for four different values of the radiation field, xuv = 10^, 
10^, 10^ and 10^ 

We have also constructed two-sided models, in which a slab of finite thickness is illu- 
minated from both sides. While the one-sided models are primarily applicable to PDRs in 
which a hot star or star cluster irradiates a nearby molecular cloud, the two-sided models 
are applicable to diffuse or translucent molecular clouds that are exposed to the ISRF. In 
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Figure 6, we present the predicted column densities as a function of the total Av,tot across 
the slab. Results are shown for clouds with density rin = 10^'^ cm~^ exposed to radiation 
with xuv = 1 and 10. In Figure 7, we present the same model predictions but with A^(H2) 
in place of Ay^tot as the independent variable. 

Figure 8 shows how the HCl and CCl"'" abundances are enhanced at elevated gas tem- 
peratures. Here we show the results of one-sided models for slabs of density nn = 10^'^ cm""^ 
exposed to radiation with xuv = 1 and 10. In addition to the standard model in which the 
temperature is determined by considerations of thermal balance, we show the abundance 
profiles for models in which the temperature is arbitrarily set to 150, 200, 300 or 400 K. 
Even temperatures as low as 150 K enhance the HCl abundance by driving the slightly en- 
dothermic reaction of CI and H2. The CCl"*" abundance is increased accordingly, a result of 
the subsequent reaction of HCl with C^. 

4. Discussion 

As expected, the results shown in Figures 3-7 indicate qualitatively-different behavior 
for the F- and Cl-bearing species. The results for F-bearing molecules are qualitatively the 
same as those obtained previously by NWS05: once H2 becomes abundant, HF becomes 
the dominant reservoir of gas-phase fluorine, with CF"*" the only other significant molecule 
containing fluorine. Because of the lower dissociative recombination rate assumed for CF"*" 
in the present work, the predicted CF+ abundances are somewhat larger. The chemistry for 
Cl-bearing molecules is somewhat more complex, by contrast, and the computer-generated 
network diagrams shown in Figure 9, together with Figures 10 - 13 in the online materials, 
help elucidate the key pathways. 

4.1. Identification of key chemical pathways 

In Figures 9 - 13, circles represent the key Cl-bearing species, with a color scale in- 
dicating their relative abundance. The color map extends from cyan to blue to purple to 
red (highest value), with grey circles denoting species that account for less than 10~^ of the 
gas-phase chlorine abundance. Arrows indicate the key reactions linking the various species: 
again, the colors of the arrows indicate the reaction rate per unit volume. The results apply 
to a one-sided model with rin = lO^cm"^ and xuv = 10^, and the different figures show 
the behavior at various depths into the cloud. The depths were chosen to illustrate several 
different regimes. 
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1) Ay = 0: At the cloud surface, there is no shielding and the abundance of molecules is 
small. The C1"^/C1 ratio is ~ 10^, being determined by the balance between photoionization 
and recombination. For the latter process, radiative recombination and charge transfer 
recombination with H - although endothermic by ~ 0.6 eV - are of roughly equal importance. 

2) Ay = 0.8: Here, H2 self-shielding has greatly increased the H2 abundance. Reaction 
with H2 is now the dominant destruction mechanism for CI"*", resulting in the formation of 
HCl"'". However, the e/H2 abundance ratio is still sufficient for dissociative recombination 
to dominate the destruction of HCl"*", although the competing reaction with H2 produces a 
small amount of H2CI"'". The dissociative recombination of HCl"'" leads to atomic CI, which 
becomes the dominant gas-phase reservoir of chlorine. Beyond Ay ~ 0.8, the abundances of 
the HCl"'" and H2CI"'' ions start to fall, along with that of C1+. 

3) Ay = 2.0: The e/H2 abundance ratio has now dropped sufficiently that HCl"*" reacts 
primarily with H2, rather than electrons, forming H2C1^. HCl is produced by dissociative 
recombination of H2CI"'", although the dominant formation process for HCl is the slightly 
endothermic reaction of H2 and CI; this reaction is still possible at the ~ 200 K temperature 
at this point within the PDR. At Ay = 2.0, CI and HCl account respectively for ~ 0.1% and 
~ 99.9% of elemental chlorine in the gas phase. 

4) Ay = 4.0: The gas temperature drops below the level at which HCl can be formed 
by reaction of H2 with CI, and the HCl abundance is lower than it was at Ay = 2.0. The 
Cl'^^/Cl ratio is now only 10~^, the photoionization rate of CI having dropped dramatically. 
Thus the formation of HCl"*" is dominated by reaction of Hj]" with CI, not H2 with Cl^. The 
Hj]" molecular ion is produced by the effects of cosmic rays; these ionize H2, forming H^, 
which then transfers a proton to H2 to form Hjj". 

5) Ay = 10.0: As at Ay = 4.0, HCl"^ is produced primarily by reaction of Hjj" with CI, 
and reaction of HCl"*" with H2 then forms H2C1"'". At this depth, the electron fraction is so low 
that dissociative recombination of H2CI"'" is no longer dominant. H2CI"'" transfers a proton to 
a neutral species of higher proton affinity than HCl (e.g. CO or H2O; this process is denoted 
by the arrow labeled "N"), yielding HCl. At this point, the external UV is almost completely 
attenuated, and the HCl destruction rate is very small. HCl now accounts for several xlO% 
of the gas-phase chlorine abundance, with atomic chlorine accounting for essentially all of the 
remainder. Our results in this regime are in good agreement with those obtained previously 
by Schilke et al (1995). H2CI"'" and CCl"*" are also present, but their abundances are just a 
few xO.01% of gas-phase elemental chlorine. HCl reacts with B.^ to form H2CI"'", although 
this rapidly leads to reformation of HCl by proton transfer to other neutral species. Several 
other destruction processes are of comparable importance for HCl destruction: cosmic ray 
induced photodissociation, and reaction with the positive ions C+ and He"*". In this regime. 
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the rates of HCl formation and destruction all scale linearly with the cosmic-ray ionization 
rate, with the result that the HCl abundance is almost independent of the latter. 

Although the diagrams shown in Figure 9-13 apply to just a single model, similar 
regimes exist for other values of xuv cind rin. Although the relevant regions move in or out 
(to larger or smaller Ay) as the assumed Xuv/f^n ratio is increased or decreased, the same 
qualitative behavior is observed. In examining similar diagrams for the entire set of models, 
we have not identified any important reaction pathway that is not apparent in Figures 9 - 
13. 



4.2. Observational implications 

To date, HF, CF"*", and HCl are the only halogen-bearing molecules to have been de- 
tected in the interstellar medium. In regard to HF, our results are identical to those presented 
in NWS05: we predict HF to be the dominant gas-phase reservoir of fluorine within both 
diffuse and dense molecular clouds; we expect the Herschel Space Observatory to detect 
widespread absorption in the HF J = 1 — transition. However, the abundances we predict 
for CF"*" lie a factor ~ 3 above the prediction of NWS05, a direct consequence of the smaller 
dissociation rate adopted for CF"*" in the present study. As noted by Neufeld et al. (2006), 
the CF"*" column densities inferred from observations of the Orion Bar were already a factor 
of ~ 4 below the predictions of NWS05, and the discrepancy between theory and observation 
is now increased to more than an order of magnitude. This disagreement may indicate that 
we have overestimated the rate coefficient for reaction of C"*" and HF. We are not aware of 
any laboratory studies of this reaction; the rate we adopt is simply the capture rate, and 
thus an upper limit on the true reaction rate. Laboratory measurements of this key reaction 
would be very desirable. 

HCl has been observed in both diffuse and dense molecular clouds. In diffuse clouds, 
ultraviolet absorption studies have led to upper limits, and in one case, a tentative detection 
toward ( Oph (Federman et al. 1995). The latter result, obtained using the Goddard High 
Resolution Spectrograph on HST, yielded an HCl column density of 2.7±1.0(1(t) x lO^^cm"^. 
Given an atomic chlorine column density of 3.0 ± 1.0(lo") x lO^^cm^^ for this sight-line 
(Federman et al. 1995), the corresponding A^(HC1)/A^(C1) ratio is 9 x 10^"^ (uncertain by 
a factor ~ 2), a value in excellent agreement with predictions of the diffuse cloud models 
presented by vDB86. Very similar results are obtained from our updated treatment of CI 
chemistry. In Figure 14, we show predictions from a series of models with xuv ranging from 
10~^ to 10^. All models apply to a slab with parameters appropriate to ( Oph, viz. density 
uh = lO^'^cm"^ and total visual extinction of 0.8 mag. The horizontal axis shows the H2 
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column density, for which the measured value is 4.2 x 10^° along this sight-line, while the 
vertical axis shows N (RCl) / N {C\) (red curve), iV(Cl+)/A^(Cl) (orange) and N(RC\+)/N{C\) 
(blue). Squares located along the curves show the results for xuv = 10~^ to 10^ from 
right to left. The H2 column density along the ( Oph sight-line requires a xuv ~ lO"'^, in 
agreement with previous studies, and the predicted 7V(HC1)/A^(C1) is in good agreement with 
the observations. The almost exact agreement between the present model and that of vDB86 
in regard to the predicted HCl abundance appears to result from the fortuitious cancellation 
of two changes in the photochemical network: the photodissociation rate for HCl and the 
photoionization rate for CI are both increased by a factor ~ 2. As in the vDB86 models, 
the predicted Cl"^ column density lies substantially below the observations. Federman et al. 
(1995) attributed this discrepancy to the presence of a significant CI"*" column density within 
HII regions along the sight-line. One important caveat must be noted. The results shown in 
Figure 14 were obtained by assuming a branching ratio of only 10% for the production of HCl 
following dissociative recombination of H2CI"''. This value was initially invoked primarily 
to explain the relative low HCl abundance derived from upper limits obtained with the 
Copernicus satellite. Models that we obtained for other branching ratios (0, 0.3 and 1.0) 
indicate that the HCl column density is linearly proportional to the branching ratio for 
values greater than 0.1 With a branching ratio of zero (i.e. with no production of HCl 
from dissociative recombination of H2C1^), a non-zero but negligible HCl column density 
(A^(HC1)/A^(C1) ~ 10~^) results from the endothermic reaction of H2 with CI. 

In dense clouds, HCl has been unequivocally detected in several sources by means of 
submillimeter observations of the J = 1—0 emission line. First detected using the Kuiper Air- 
borne Observatory toward Orion (Blake, Keene & Phillips 1985) and then Sgr B2 (Zmuidz- 
inas et al. 1995), the J = 1 — transition has also been observed using the ground-based 
Caltech Submillimeter Observatory (CSO) under good atmospheric conditions. CSO obser- 
vations have led to additional detections toward Orion A, Mon R2 (Salez, Frerking & Langer 
1996) and several other sources (Phillips et al. 2009), as well as mapping observations to- 
ward OMC-1 (Schilke, Phillips & Wang 1995). Typical column densities derived for these 
dense clouds lie in the few xlO^'^ to few xlO^^cm"^ range, in reasonable agreement with the 
predictions of our model for regions at high density exposed to strong UV radiation (Fig. 
5). As noted in the observational studies cited above, depletion plays an important role in 
limiting the fractional abundance of HCl. 

Our study identifies two additional Cl-bearing species that are potentially detectable: 
the molecular ions HCl"*" and H2CI"'". These ions are isoelectronic with OH and H2S re- 
spectively. They are most abundant near cloud surfaces, where the photoionization rate 
for HCl is highest, and show column densities that are an increasing function of Xuv/nn- 
Thus, PDRs subject to strong UV irradiation present attractive targets for searches for 
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HCl"'" and H2C1^ at millimeter and submillimeter wavelengths. For example, in a PDR with 
TT-H = lO^cm"'^ and xuv = 10^, column densities of 5.7 and 2.6 x 10^^ cm~^ are predicted for 
HC1+ and 13201"^. These values apply to sight-lines perpendicular to the illuminated surface; 
in edge-on PDRs, they can be significantly enhanced by limb-brightening. The rotational 
spectrum of HC1+ has been measured using laser magnetic spectroscopy (Lubic et al. 1989). 
The lowest-lying rotational transition, the ^n3/2</ = 5/2 3/2 multiplet near 207.6/im, lies 
in a wavelength range accessible to the HIFI instrument on the Herschel Space Observatory. 
The rotational spectrum of H2CI+ has also been the subject of a laboratory study at high 
spectral resolution; here, the measurements of Araki et al. (2001) indicate that the lowest- 
lying rotational line of ortho-H2Cl^, the lio — loi transition, lie near 189.2 and 188.4 GHz 
respectively for the Hg^Cl"'' and H2^C1"'' isotopologues. (These transitions are split into 6 
and 4 hyperfine components, covering ~ 56 and 14 MHz respectively.) This spectral region 
lies in the wing of the strong 183 GHz telluric water absorption line, but is nevertheless 
observable from ground-based observatories under favorable - but by no means unusual - 
atmospheric conditions. In one of the earliest studies of interstellar chlorine chemistry, the 
possibihty of detecting HC1+ in diffuse molecular clouds was discussed by Jura (1974), who 
considered the A^E^ ^ X^H band that is accessible in the near-ultraviolet region. Given 
an oscillator strength of 4.15 x 10""^ for the strongest vibrational band (Pradhan, Kirby & 
Dalgarno 1991), and given the abundances predicted in Figure 17, our results for C, Oph 
would imply equivalent widths of only ~ 4 x 10~^ A, well below the limit of detectability. 

Finally, we note that, unlike CF"*", the CC1+ ion is typically predicted to have an ex- 
tremely low abundance. Although CCl"'" is produced rapidly by reaction of HCl and C"*", 
these two species show very little overlap: the HCl abundance is small in the surface layers 
where C"*" is abundant, unless the temperature is increased by shock heating, and the C"^ 
abundance is small in the deep cloud interiors where HCl is relatively abundant. Enhanced 
CCl"*" abundances are possible when shock heating is present along with UV irradiation. 

We are very grateful to W. D. Geppert and M. Hamberg for informing us about recent 
laboratory results on H2CI+ dissociative recombination, in advance of their publication. We 
thank Paule Sonnentrucker and Paul Bryans for several useful discussions. We gratefully 
acknowledge the support of a grant from the NASA Herschel Science Center's Theoretical 
Research/Laboratory Astrophysics Program. 
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Table 1. Reaction list for Cl-bearing species 



Reaction Rate or rate coefficient 



CI + Ho ^ HCl + H 



HCl + H2 ^ CI + H2 



C+ + HCl - 
H+ + HCl - 
He+ + HCl 
C1+ + H2 - 
HC1+ + H2 
CC1+ + e - 
HC1+ + e - 
H2CI+ 
H2CI+ 
H2CI+ 
H2CI+ 
Cl + H 



CC1+ + H 
H2 + H2CI+ 
H + C1+ + He 
HC1+ + H 

H2CI+ + H 
C + Cl 
H + Cl 
^ HCl + H 
> CI + products 
^ HCl + HCO+ 



e — 
e — 
CO 

H2O ^ HCl + H3O- 
H2CI+ + H 



HCl + hi^ H + CI 



HCl + hu HCr 



exp(-2214K/r) cm^ s^^ 
12 ( r/300 K)i-^2 (_i544 K/ T) cm^ s'^ 



1.49 X 10-^^ exp(-1763K/r) cm^ s-^ 



2.52 X 10"" 

2.86 X 10 

-1-11 

1.69 X 10-12 ( r/300 K)i-72 

1.84 X 10"^ (T/300K) 
6.35 X 10-9 (T/300K)-o-202 

5.51 X 10-9 (T/300K) 
1.0 X lO-^cm^s-i 
1.3 X lO-^cm^s-i 



exp(-1093K/r)cm3 s^ 
1 



(T < 354 K) 
(T > 354 K) 
(T < 354 K) 
(r > 354 K) 



0.202 ^^3 g 



cni"^ s 1 
0-202 cm^s-i 



cm s 

n ^ - 
cm s 



cm'^ s 1 



—1 

cm s 



0.194 g 



2.0 X 10-^ (r/300 K) 
2.0 X 10-^(T/300K) 
1.2 X 10-s (r/300 K)-o-S5 
1.08 X 10"^ (r/300 K)-o-85 
7.8 X 10-1° cm^ s-i 
3.7 X 10-9 (r/300 K) 
1.0 X lO-^cm^s-i 

1.7 X 10-9 xuv[lE2{2.1Av) + lE2{2.1[Av,tot ~ Ay])] s'^ 

+ 1370 Cp /(I -^) 
1.0 X 10-10 xuvl i^2(3.45Av) + 1^2 (3.45 [v4y,tot - Ay])] s"! 

X exp(-iV(H2)/3.5 x 102icm-2)/(l + [iV(H2)/3.2 x 1020cm-2]V2) 

+ 185Cp/(l-a;) 



Table 1 — Continued 



Reaction Rate or rate coefficient 



CI + hi/ ^ C1+ + e 4.8 X 10-" xuv[ |^2(3.6Ay) + 1^2(3.6 [A v,tot - Ay])] 

X exp(-A^(H2)/4.28 x 10^^cm-^)/{l + [N{R2)/4:.9 x lO^Ocm-^] V^) 
+ QlCp/{l-uj) 

C1+ + e ^ CI + hu 1.34 x 10^" (T/300 K)"^-'^^^ cm^ s^^ 

C1+ + H ^ C1 + H+ 6.2 X 10-" (T/300 K)°-^9 exp (-6920 K/T) cm^ s-^ 

CI + H+ ^ C1+ + H 9.3 X 10-" (T/300 K)^-^^ (_232 K/ T) cm^ s-^ 

HCl + H+ ^ HC1+ + H 3.3 X 10-^ (T/300 K)i °o cm^ s-^ 
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Table 2. Reaction list for F-bearing species 



Reaction 



Rate or rate coefficient 



F + H2 

F + H2O 

C+ + HF - 
Si+ + HF - 
H+ + HF - 
He+ + HF 
F+ + H2 - 

HF+ + H2 
CF+ + e - 
SiF+ + e - 
HF+ + e - 
H2F+ + e - 
H2F+ + e - 
F + CH^ 
F + OH^ 
F + H+- 
HF + hi/ 



HF + H 

HF + O*^ 

CF+ + H 
SiF+ + H 
H2 + H2F+ 
^ H + F+ + He 
F 



H 



^ HF+ 
^ H2F+ + H 
C + F 
Si + F 
^ H + F 

HF + H 
-> F + products 
HF + C 
HF + 
H2F+ + H 
H + F 



s —1 
cm s 



0.15pj^3g 1 



1.0 X 10-1° [exp (-450 K/T) + 0.078 exp (-80 K/T) 

+ 0.0155 exp (-10 K/ T)] cm^ s-^ 

1.6 X 10-11 cm^s-i 

7.2 X 10-9 (T/300K)-° i5cm3s-i 

5.7 X 10-9 (T/300 K)-°-i5 cm^ s'^ 

1.2 X 10-® (T/300 K)-°-i5 

1.1 X lO-s (T/300 K) 
6.24 X 10-1° cm^' s-i 

3.8 X 10-i°cm^s-i 

1.3 X lO-^cm^s-i 

5.2 X 10-* (T/300 K)-° « 

2.0 X 10-^ (T/300 K)-°-5cm3s-i 

0.5 g-i 

s —1 
cm s 

0.5 g-i 



cm s 



2.0 X 10-^ (T/300 K) 
3.5 X 10-^ (T/300 K)-°-5 

3.5 X 10-^ (T/300 K) 

1.6 X lO-i^cm^s-i 
1.6 X 10-i°cm3s-i 
4.8 X lO-i^cm^s-i 

1.17 X 10-10 Xf/y[|^2(2.21Ay^ 
+ 124Cp/(l-cu) 



lE2{2.21[Av,tot-Av])]s-' 



^Bold face denotes those reactions for which the adopted rates are different from those assumed 
by NWS05 
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H2 

4.48 
















LiH 
2.41 


BeH 
2.24 


BH 
3.44 


CH 
3.49 


NH 

3.22 


OH 
4.39 


HF 
5.87 




NaH 
2.04 


MgH 
1.99 


AIH 
2.95 


SiH 
2.98 


PH 
2.87 


SH 
3.65 


HCI 
4.43 




KH 
1.87 


CaH 
1.77 


GaH 
2.78 


GeH 
2.73 


AsH 
2.66 


SeH 
3.11 


HBr 
3.76 





ScH 


TIH 


VH 


CrH 


MnH 


FeH 


CoH 


NIH 


CuH 


ZnH 


2.09 


2.08 


1.67 


2.27 


1.27 


1.59 


1.97 


2.57 


2.64 


0.83 



Fig. 1 - Dissociation energies, Dq (at K in eV), for the neutral diatomic hydrides (derived 
from values recommended in the NIST Chemistry Web book and in Armentrout & Sunderlin 
1992, Chen, Clemmer & Armentrout 1991, 1993; Binning & Curtiss 1990; and Huber & 
Herzberg 1979) 
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Fig. 2 - Dissociation energies of HX and HX+ for X = H, C, N, O, F, Si, P, S, and CI. 
Atoms with IP < IP(H) appear in blue; those with IP > IP(II) appear in red. The dashed 
lines indicate the dissociation energy of H2 (4.48eV). 
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Av (mag) 

Fig. 3 - Abundance profiles for a one-sided slab model with nn = 10^ cm~'^ and xuv = 10^. 
Solid lines show the results for an assumed primary cosmic ray ionization rate of 
1.8 X 10"^'^ s~^ per H nucleus (standard case). Dashed lines show the results for an assumed 
primary cosmic ray ionization rate of 1.8 x 10~^^s~^ per H nucleus. 
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Fig. 4 - Profiles of temperature and F- and Cl-depletion factor, for a one-sided slab model 
with nn = 10^ cm~'^ and xuv = 10^. Solid lines show the results for an assumed primary 
cosmic ray ionization rate of 1.8 x 10^^'' s^^ per H nucleus (standard case). Dashed lines 
show the results for an assumed primary cosmic ray ionization rate of 1.8 x 10~^^ s~^ per H 
nucleus. 
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5 - Molecular column densities, for one-sided slab models 
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Fig. 6 - Molecular column densities, for two-sided slab models, as a function of Ay across 
the slab 
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Fig. 7 - Molecular column densities, for two-sided slab models, as a function of A^(H2) 
across the slab 



- 29 - 




Fig. 8 - Effect of enhanced temperatures on the abundances of HCl (red curves) and CCl"'' 
(green curves). Results of one-sided models for slabs of density tt-h = 10^'^ cm~^ exposed to 
radiation with xuv = 1 (left panel) and 10 (right panel). In addition to the standard model 
in which the temperature is determined by considerations of thermal balance (solid curve in 
each set), we show the abundance profiles for models in which the temperature is arbitrarily 
set to 150 K (dotted), 200 K (dashed) , 300 K (dot-dash), and 400 K (triple-dot-dash) 
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CI abundance = 1 .8c — 7 / K/cx. rcccticn rate = 8.5c— 14 cnn '' 



-8 -6 -4 -2 

bg(reaction rate/rTox. rate) or log(Qbundance/CI abundance) 

Fig. 9 - Chemical network diagram, for a one-sided slab model with hh 
Xuv = 10"^. The figure is labeled by the depth, Ay, in visual magnitude. 



10^ cm ^ and 
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CI abundance = 1 .8c — 7 / K/ax. rcccticn rate = 4.3c— 12 cnn 
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bg(reaction rate/rTox. rate) or log(Qbundance/CI abundance) 



Fig. 10 - Chemical network diagram, for a one-sided slab model with nn = 10^ cm ^ and 
Xuv = 10"^. The figure is labeled by the depth, Ay, in visual magnitude. 





CI abundance = 1 .8c — 7 / Wax. rcccticn rate = 5.6c— 14 cnn 
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bg(reaction rate/rTox. rate) or log(Qbundance/CI abundance) 



Fig. 11 - Chemical network diagram, for a one-sided slab model with nn = 10^ cm ^ and 
Xuv = 10"^. The figure is labeled by the depth, Ay, in visual magnitude. 





CI abundance = 1 .8c — 7 / Wax. rcccticn rate = 2.9e— 17 cnn 
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bg(reaction rate/rTox. rate) or log(Qbundance/CI abundance) 



Fig. 12 - Chemical network diagram, for a one-sided slab model with nn = 10^ cm ^ and 
Xuv = 10"^. The figure is labeled by the depth, Ay, in visual magnitude. 
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A,= 1 0.00 




CI abundance = 1 .8c — 7 / K/ax. rcccticn rate = 5.8c— 16 cnn '' 
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Fig. 13 - Chemical network diagram, for a one-sided slab model with nn 
Xuv = 10"^. The figure is labeled by the depth, Ay, in visual magnitude. 
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Fig. 14 - Predictions from a series of models with xuv ranging from 10~^ to 10^. All models 
apply to a slab with parameters appropriate to ^ Oph, viz. density uh = 10^'^ cm~^ and 
total visual extinction of 0.8 mag. The horizontal axis shows the H2 column density, for 
which the measured value is 4.2 x 10^° along this sight-line, while the vertical axis shows 
iV(HCl)/iV(Cl) (red curve), A^(Cl+)/iV(Cl) (orange) and iV(HCl+)/iV(Cl) (blue). Squares 
located along the curves show the results for xuv = 10~\ 1, 10 and 10^ from right to left. 



